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Abstract. In this paper, periodic homogenization of a steady fluid flow in 
fissured porous solids with imperfect interfacial contact is performed via two- 
scale asymptotic method. 



1. Introduction 

Micro-models for the fluid flow in fissured porous medium made of two inter- 
acting porous systems with different permeabilities are very important in reservoir 
petroleum, civil engineering, geophysics and many other areas of engineering. In 
this paper, we shall deal with the homogenization of fluid flow in fissured porous 
media. There are many papers devoted to the subject, see for instance [H El [7] 
and the references therein. Here, we shall be concerned with periodic porous media 
made of two components with highly contrast of permeabilities and an imperfect 
contact on the interface of these solids. One of the two porous structures is asso- 
ciated with the fissures and the other one with the porous block or matrix. More 
precisely, we shall consider a Newtonian fluid flow, in a periodic fissured porous 
medium ft = fif U $71 U F*^, where fif is the fissures region, $71 is the block and 
r*^ is the interface that separates these two regions. The micro-model is based on 
mass conservation for the fluid in each phase, combined with the Darcy's law . We 
assume that fif and r2| are in imperfect contact. That is, an exchange flow barrier 
formulation on the interface will be considered, see for instance, [TJIH1I2]- The 
macro-model is derived by means of two-scale asymptotic technique {Cf. HH]) 
and jusitifled with the help of the two-scale convergence method [3] • It will be seen 
that the macro-model is, in some sense, the limit of a family of periodic micro- 
models in which the size of the periodicity approaches zero. It is shown that the 
overall behavior of fluid flow in such micro-scale heterogeneous media, is a stan- 
dard one. Namely, the macroscopic equation is of elliptic type and it has the form 
— div(A''Vu) = F. The only novelty here is that the fluid flow presents an ex- 
tra source surface density of the exterior boundary. We mention that no double 
porosity effects occur at the macroscale description. 

The paper is organized as follows: Section 2 is devoted to the problem setting 
of the micro-model. In Section [31 we shall be concerned with the derivation of the 
homogenized model via the formal procedure of two-scale asymptotic ananlysis. 
Section |4] is devoted to the rigorous derivation of the homogenized model obtained 
in the previous section by the two scale convergence method. Finally, we end this 
paper with some comments and conclusion. 
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2. Setting of the Problem and the main result 

We consider il a bounded, connected and smooth open subset of {N > 2) . 
We assume that fl is made of a large number of identical cells eY, where e > is 
a sufRciently small parameter (e ^ 1) and Y —]0, denotes the generic cell of 
periodicity. Let Y be divided as follows: Y := Yi UF where Yi, Y2 are two sub- 
domains of Y and F is the interface between Yi and Y2 . Namely F = Yi n I2 ■ We 
denote v the unit normal of F, outward to Yi. Let xi and X2 denote respectively 
the characteristic function of Yi and Y2, extended by Y-periodicity to M.^ . For 
i = 1 or 2 we set 

ni ■.^{xen■.x^{^) = l}. 

Let F"^ Hf n Hf . We will assume that II| C f^, so that = F"^ and d^l = 
dfl U F*^. Let Zi — Uj,gz" O^i + k). As in [3^, we shall assume that Zi is smooth 
and connected open subset of . 

Let A and B denote respectively the permeability of the medium Zi and Z2. 
We assume that A (resp. B) is continuous on M^, Y— periodic and satisfies the 
ellipticity condition: 

AC • C > C|er, (resp. i?e ■ e > C e 

where, here and in what follows, C denotes various positive constants independent 
of £. Let fi be a measurable function representing the internal source density of 
the fluid flow in Vlf and let h be the non-rescaled hydraulic permeability of the thin 
layer F*^. We suppose that fi e (fi) and that h is a. continous function on , 
Y— periodic and bounded from below: 

h{y)>ho>0 X E R^. 

To deal with periodic homogenization with microstructures, we shall denote for 
X G M^, A^ (x) = v4 (f ), B^ {x) = B (f ) and (x) = eh (f). 

The micro-model that we shall study is given by the following set of equations: 



(2.1a) - div {A'Vu') = fi in 17f , 

(2.1b) - e^div {B^Vv") = /a in rj|, 

(2.1c) A^Vu" ■ = e^B'S/v" ■ on F^ 

(2. Id) A^Vu" ■ = -eh" (u^ - v^) on F^ 
(2.1e) = on dVl 



where stands for the unit normal of F*^ outward to Of. It is obtained by the 
y— periodicity extension of v. Let us mention that fif represents the fissured space 
region with permeability A and f^l the block region with permeability e^B. The 
quantitities u*^ and are the fiuid flow velocities in fJf and respectively. Note 
tha we have chosen a particular scaling of the permeability coefficients in (|2.1bp . 
This means that both terms \'SIu'^\ dx and \Vv'^\ dx have the same order 

of magnitude and thus leading to a balance in potential energies. Equations (|2.1ap 
and (|2.1bl) express the conservation of mass of fluid with Darcy's law in fif and 
respectively. The first equation describes the flow in the fissured regions with 
large permeability and the second describes the flow in the block system regions 
with low permeability. For more details, we refer the reader to Arbogast, Douglas, 
and Hornung [5] (see also Allaire [3]). Condition (I2.1cp expresses flux continuity 
across F^, the interfacial flow thin layer with permeability given by . Condition 
(|2.1d|) models the imperfect contact between the block and the fissures [HI [H] ■ The 
condition (|2.1ep is the standard homogeneous Dirichlet condition on the exterior 
boundary of fi. 
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To set the mathematical framework of om- Porblem (|2.1ap - (|2.1ep . we need to 
introduce the following space = {H^ {0.1) n (fi)) x The space 

is equipped with the norm 



||(^, = W^vWl^i^m) + ^ II V^llL2(n|) + £ 11^ - V-lll^tre) . 

The weak formulation of (j2.1a[) - (j2.1ep can be read as follows: find {u'^,v^) G 
such that for all v = {ip, ip) G , we have 



(2.2) / {v^ - v^) {(f - V) ds"" {x) = [ fupdx + I f2^dx. 



where dx and ds"^ {x) denote respectively the Lebesgue measure on and the 
HausdorfF measure on F^. 

In view of the assumptions made on Ai, fi and h, we can easily establish the 
following existence and uniqueness result whose proof is a slight modification of 
that given by H. I. Ene and D. Polisevski [5] and therefore will be omitted. 

Theorem 2.1. For any sufficiently small e > 0, there exists a unique couple 
{u^,v^) £ , solution of the weak problem such that 

(2.3) Uu^,vn\\H^ <C. 

Now, thanks to the a priori estimates (|2.3p . one is led to study the limiting 
behavior of the sequence {u'^,v'^) as e approaches 0. This is summarized in the 
main result of the paper: 

Theorem 2.2. Let {u'^,v'^) G H'^ be the solution of the weak system 112. S^) . Then, 
up to a subsequence, there exists u G Hq (f2) and vq G (^O; (^^2)^ such that 

u'^ u in I? ip.) weakly, 

Vq {y) dy in L'^ {n) weakly. 

Let = xi (f ) ■'^'^ + X2 (f ) v'^ denote the overall pressure. Then, w'^ weakly con- 
verges to U — u + G where U is the unique solution to the homogenized model: 

-div (A^VU) =F m O, 

u^G on an. 

Here, A^ (resp. F, G) is given by JOT]) (resp. |0^ , ^0^) ). 

To prove this theorem, we shall first use the two-scale asymptotic procedure to 
formally derive the homogenized model and apply then the two-scale convergence 
technique to rigorously justify the homogenization. This is the scope of the next 
two Sections. 



3. The formal homogenization procedure 

The purpose of this section is to formally construct the homogenized system of 
(l2.1ap - (l2.1eL via the two-scale asymptotic expansion method ^ ilOj. Assume the 
following ansatz: 

(3.1) (x) = UQ{x,y) + eui{x,y) + e'^U2ix,y) + . . . , 

(3.2) (x) = vo{x,y) + svi{x,y) + s^V2ix,y) + . . . 
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where y — x/e, the unknowns uq, vq, ui,vi, U2,V2.-- are y-periodic in the second 
variable y. Plugging the above expansions p.ip - p.2p into the set of equations (j2.1al - 
I2.1ep and identifying powers of e yields a hierarchy of boundary value problems. 
At the first step, it is not difhcult to observe that Equation (|2.1bl) at £~^, and 
Equation (|2.1cp at e^^ orders yield that uq is independent of y e Fi. That is 

Mo (a;, y) = u (x) . 

Next, Equation (j2.1bp at and Equation (|2.1cp at e° orders show that the cor- 
rector term ui may be written as 




where, for 1 < j < iV, Ljj e (Yi) /R is the unique solution to the following cell 
problem: 

/ {AVyLJ„VyC)dy^ f {-Aej,VyC)dy, C^H'#iY,), 

JYi JYi 

where (e^) is the canonical basis of M^. In p.3p . u{x) is any additive constant. 
At the final step. Equation ((2?Ta| - (j2.1bp at e° and Equations ((2?Tc)) - (j2.1dp ) at 
orders yield the system: 

(3.4) - divj^ {A\7yU2) = fi + divy {AV^ui) + div^ {A (V^ui + V^u)) in 17 x Yi, 



(3.5) - divj, {BVyVo) = /2 in 17 X Y2; 

(3.6) AVyU2 ■ V = —AVxUi ■ V + BVyVo ■ ly on fl xT, 

(3.7) A\7yU2 ■ V = — AVajMi ■ V — h{u — wq) on 12 x F, 

(3.8) y^vo (x, y) , U2 [x, y) Y - periodic. 



The weak formulation of Equations ((Xi)) . ((X71) - (P^ is 

/ [A\IyU2.VyOdy^{F,C). Cei?#(n), 

JVi 

where 

{FX) = / div,(yl(VyMi + Vu))C 

JYi 

- [ ^V,uiVyC+ / h{u-vo)Cds{y)+ [ M 
JYi Jr JYi 

where ds (y) denotes the Hausdorff measure on F. Again, using the Divergence The- 
orem ( as in lOJ ) , a necessary and sufficient condition for the existence and unique- 
ness of U2 G (Yi) /R is that F satisfies the compatibility condition {F, 1) = 0. 
This reads in our case as follows: 

(3.9) / {divx{A{VyUi + Wu)))dy+ I h{u-vo)ds{y) = I fidy. 
JYi Jr Jyi 

Using (|3.3p . Equation (|3.9p becomes 

(3.10) -div(^'^Vu)+ [ h {u - vo) ds (y) = [ fidy 

Jr JYi 

where A^ = ("y)i<ij<Ar given by 

(3.11) = / ^ i^y^^ + ' C^y^j + ej) dy. 

JYi 
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On the other hand, it is easy to see that vq can be written as 

(3.12) voix,y)-u{x) = a{y)f2ix), (x, y) G 17 x Fa 
where a G {Y2) is the unique solution of the following problem 

r -divy {BVya) = 1 in ^2, 

(3.13) <^ 

I B\7ya ■ V — ha = on F. 

Therefore, Equation p.lOp yields 

(3.14) - div (A^Vu) = \Yi\fi + a/2 F* 

where d = Jp ads (y). The boundary condition for u is obtained from (j2.1ep at 
order and it reads 

(3.15) u = on dn. 

Let us observe that the overall pressure = xi (-) u^ + X2 (-) two scale con- 
verges to u + X2Ctf2 (see Definition 14. 1 1 below) . Consequently, w'^ weakly converges 
to U = u + G where 



(3.16) G=(^y^aj/2. 

Equation (j3.14p is the so-called macroscopic equation for u. Moreover, if G is 
sufhciently smooth, then the homogenized model for the weak limit U is as follows: 

(3.17) -div (A^VC/) = F in 17, 

(3.18) C/ = G on dVL 

where F = F* +div (A'^VG) . Note that in the homogenized problem (PT7)) - (|XT5)) . 
a non homogeneous Dirichlet is derived, on the contrary to the micro-model, where 
homogeneous boundary condition is prescribed, see condition (j2.1ep . This extra 
source surface density essentially arises from the fact that 

(1) Blocks have low permeability; 

(2) non null source volumetric density on the blocks; 

(3) and the contribution of the fluid flow in the vicinity of the thin layer where 
the contact is assumed imperfect. 

4. Two-scale convergence approach 

In this section, we will derive the homogenized system of (j2.1ap - (j2.1ep by the 
two scale convergence method, see j2|. First, we define C^iY) to be the space 
of all continuous functions on which are F-periodic. Let the space {Y) 
(resp. (Fi), i = 1, 2) to be all functions belonging to Lf^^ (R^) (resp. L^^^ (Zi)) 
which are F-periodic, and (Y) (resp. {Yi)) to be the space of those functions 
together with their derivatives belonging to L"^ (Y) (resp. (^i))- Now, we recall 
the definition and main results concerning the method of two-scale convergence. For 
full details, we refer the reader to SIH]. 

Definition 4.1. A sequence (w^) in (17) two-scale converges to v E {H. x Y) 
(we write "-^^ v) if, for any admissible test function Lp E (Jl;C^{Y)) , 



lim 



1 / v'^ {x) Lp (x . —\ dx ^ / V {x , y) If {x , y) dxdy 
In V e/ InxY 



Theorem 4.1. Let (v^) be a sequence of functions in L^(S7) which is uniformly 
bounded. Then, there exist v G {ft x Y) and a subsequence of (w^) which two- 
scale converges to v. 
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Theorem 4.2. Let (v^) be a uniformly bounded sequence in H^{fl) (resp. HQ{n)). 
Then there exist v G (ft) (resp. Hq{Q)) and v G L'^{il; H^{Y)/R) such that, up 
to a subsequence, 

v; Vv" Vv + Vyi). 
Theorem 4.3. Let (w^) be a sequence of functions in H^(fl) such that 

\\v'\\mn)+^\\^^"\\mnf<C- 

Then, there exist v G L (n:H^{Y)^ and a subsequence of (f^), still denoted by 
(w^) such that 

and for every if (f2; C#(y)), we have: 



e->0 



lim / ev^ 



{x) LP (x, - ) ds" [x) = I V (x, y) (f {x, y) dxds [y) . 



Here and in the sequel ds (y) denotes the Hausdorff measure on T. 

As a direct application of the theorems listed above (Thms 14. 1114. 3p and the a 
priori estimates (|2.3L we give without proof the following two-scale convergence 
result concerning the solutions {u^ , u^) of the Problem (12.21) . 

Theorem 4.4. There exists a subsequence of {u"^ , v^), solution of \2/0\) . still denoted 
and there exist unique u G H}^ (fi), ui G L^{fl;H^ (Y) /R) and vq G 



such that 



(4.1) Xiu Xiu, X2V X2V0 



an 



d 



(4.2) xfVu^ XI {^U + VyUi) , £X2Vf" X2yyVo- 
Moreover, the following convergence holds: 

(4.3) lim / e{u' -v')ip'ds^x) ^ [u - vq) i^dxds {y) 



e-i-O 



for any ijj G 2? (fJ; {Y)) with tp^ (x) — (x,x/e). 

To determine the limiting equations of the system (|2.2p . we begin by choosing 
the adequate admissible test functions. Let ip'^{x) = ip {x)+eipi (^x, — ^ and 'ip'^{x) = 

■0 (x, where (p e V (fl) and ipi, ip £ V (d,; (F)^ . Taking (/? = and -0 = "0" 



£ 

in (12. 2L we obtain 



/ Ai^-^Vu" {\/(p + \/y(pi)dx+ eB (^-^Vv'^'Vytpdx + 
(4.4) ej h{^^{u' -v'){ip-i))ds'' {x)+eR^ ^ j fupdx+j f2ipdx 



where 



= J ^(^) Vu^V:,(^i (x, J) dx + e^ B V«^V^V (a;, J) dx 
ef /i(^)(ii^-t;^)^i(x,J)ds=(x). 
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In view of (|4.2p . one can deduce that 

lim / ^{~) (^'^ (^) + ^a"^! f) ) '^^ = 

Xi (y) A (y) (Vu + Vj,ui) {\7(p (x) + Vy(pi {x, y)) dxdy 



e->-0 



nxY 



and 



-0 4'''©^"'^^^ = 
X2 (y) B (y) Vvo^Jyi^ (x, y) dsdy. 

By virtue of (|4.3p . we find that 

lim £ y 1^-^ (m' - i;"^) (a;) - ijj [x, ^ ds^ (x) = 

/i (y) (u - Wo) (95 {x) - -0 (a;, y)) dxds (y) . 

Oxr 

We observe that = O (1) and, by cohecting all the preceeding limits, we get the 
following limiting equation of (|2.2I) : 

(y) (Vu + VyUi) (V(/3 + Vy<y9i) dxrfy + / B {y)V yVQ\/ yipdxdy + 
rtxYi JnxY2 

h (y) {u - vq) {ip - -0) dxds (y) 

sixr 

(4=5) / hipdx+ / /20da;. 
JnxYi JnxY2 

By density argument, the equation (|4.5p still holds true for any ((/?, (pi, (/?2) G 
(H) X £2 l^f^; iJi, (^1) /k) X L2 (^f^; iJ^, (^2)^ We can summarize the preceding 
by observing that these equations are a weak formulation associated to the two-scale 
homogenized system (|4.6p - (|4.12p below: 

(4.6) -divy {A (V-u + V^ui)) = a.e. inflxYi, 

(4.7) -div J, (i?Vj;Uo) = /2 a.e. in Q XY2, 

-div A{Vu + '\/yUi)dy 

(4.8) J h (y) [u - vo] ds (y) = /i a.e. in n, 

with the transmission and boundary conditions: 

(4.9) {A (Vu + VyUi)) • 1^= a.e. on x r, 

(4.10) BVyVQ ■ v= - h (y) [u - vq] a.e. on il x F, 

(4.11) yi — 5- ui, uo i^-periodic, 

(4.12) u = on dn. 

Let us first note that Equations (|4.6p and (|4.9p lead to the relation (j3.3p . Sim- 
ilarly, Equations (|4?7l) . (|4TT0)) and (|4T1|) yield (l3TT2|) . Equation (|48)) is the same 
as (|3.9p . With these remarks in mind, one can easily recover the homogenization 
procedure of the preceeding section and thus we have rigorously justified the formal 
two-scale asymptotic expansion method. 
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5. Conclusion 

We have used the homogenization theory to derive a macro-niodel for fluid flow 
in fissured media with microstructures, in which blocks present very low perme- 
abilities. Moreover, the contact between the block region and the fissures region 
is of imperfect type. The micro-model that we considered in this paper is with 
a homogeneous Dirichlet boundary condition prescribed on the exterior boundary. 
We have shown that the overall behavior of fluid flow in such heterogeneous media 
with low permeability at the micro-scale is a classical problem except that in the 
vicinity of the surface, there is an additional source density arising from the source 
density of the blocks at the micro-scale and the fact that the contact at micro-scale 
is imperfect. 
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